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^ C*| This paper presents the Euler-Lagrange equations for fractional variational problems with 



multiple integrals. The fractional Noether-type theorem for conservative and nonconservative 
generalized physical systems is proved. Our approach uses well-known notion of the Riemann- 
Liouville fractional derivative. 
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1 Introduction 



Fractional Calculus is the branch of mathematical analysis that studies derivatives of arbitrary 
(real or complex) order. Fractional differentiation and fractional integration are now recognized 
as vital mathematical tools to model the behavior and to understand complex systems, classical 
or quantum, conservative or nonconservative, with or without constraints. Several books on the 
subject have been written, illustrating the usefulness of the theory in widespread fields of science 
and engineering (see [35] for a review). Some researchers refer to FC as the calculus of the XXI 
century. The fractional operators are nonlocal, therefore they are suitable for constructing models 
possessing memory effect. 

The fractional calculus of variations unifies the calculus of variations and the fractional calcu- 
lus, by inserting fractional derivatives into the variational functionals. This occurs naturally in 
many problems of physics, mechanics, and engineering, in order to provide more accurate mod- 
els of physical phenomena. The physical reasons for the appearance of fractional equations are, 
in general, long-range dissipation and nonconservatism. A theory of the fractional calculus of 
variations started in 1996 with works of Riewe |35) . in order to better describe nonconservative 
systems in mechanics. He explained: "For conservative systems, variational methods are equiv- 
alent to the original used by Newton. However, while Newton's equations allow nonconservative 
forces, the later techniques of Lagrangian and Hamiltonian mechanics have no direct way to deal- 
ing with them." Riewe formulated a Lagrangian with fractional derivatives, and obtained the 
respective Euler-Lagrange equation, combining both conservative and nonconservative cases. Re- 
cently, several approaches have been developed to generalize the least action principle and the 
Euler-Lagrange equations to include fractional derivatives. Investigations cover problems depend- 
ing on the Riemann-Liouville fractional derivative (see, e.g., [5J HH I3Z]), the Caputo fractional 



derivative (see, e.g., [H 1251 [52]). and others [Ml H3 GIB] ■ We cite here also [TO], where an in- 
teresting approach known as the Fractional Action-Like Variational Approach (FALVA) to model 
nonconservative dynamical systems is presented, and its extension counterparts [131 1141 II 51 116] . 
In (TTJ [T2] the authors study multidimensional FALVA problems, which are generalization of 1D- 
FALVA problems with fractional derivatives defined in the sense of Cresson. We also mention 
study of multidimensional fractional variational problems in context of the coherence problem 
7, 8 or via the fractional derivative defined in the sense of Jumarie [2]. In the present manuscript 
we rather generalize the Riewe approach to multidimensional fractional variational problems and 
derive corresponding Euler-Lagrange equations. 

Conservative physical systems imply frictionless motion and are a simplification of the real 
dynamical world. Almost all processes observed in the physical world are nonconservative. For 
nonconservative dynamical systems the conservations law are broken so that the standard La- 
grangian or Hamiltonian formalism is no longer valid for describing the behavior of the system. 
However, it is still possible to obtain the validity of Noether's principle using the fractional calcu- 
lus of variations. Roughly speaking, one can prove that Noether's conservation laws are still valid 
if a new term, involving the nonconservative forces, is added to the standard constants of motion. 
The seminal work p2] makes use of the notion of Euler-Lagrange fractional extremal introduced 
by Riewe, to prove a Noether-type theorem that combine conservative and nonconservative cases. 
After that, other results about the subject appeared [5J [71 [TSJ [TS]. Multidimensional version of 
Noether's theorem is of most relevance to modern physics, particularly to relativistic field theories 
and to theories of gravity. To the best of our knowledge, the fractional Noether-type theorem 
for multidimensional Lagrangians is not available. Such a generalization is the aim of this paper. 
We prove the fractional Noether theorem for multiple dimensional Lagrangians, and we show how 
our results can be use. We trust that this paper will open several new directions of research and 
applications. 

The paper is organized as follows. In Section [2] we review the necessary notions of fractional 
calculus. Our results are given in next sections: in Section [3] we derive the Euler-Lagrange 
equations for fractional variational problems with multiple integrals (Theorem[T]), these conditions 
are then applied to physical problems; in Section [4] we prove the fractional Noether-type theorem 
(Theorem (2). 

2 Fractional Calculus 

In this section we review the necessary definitions and facts from the fractional calculus. For more 
on the subject we refer the reader to [23" 1 134 ) . 

Let [a, b] be a finite interval, a 6 R and < a < 1. We begin with Riemann-Liouville Fractional 
Integrals (RLFI) of order a of function / which are defined by: the left RLFI 



J°f(x) = -i- [ X (x - tr~ l f(t)dt, x>a, (1) 

1 (a) Ja 



the right RLFI 
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where F(-) represents the Gamma function. For a = 0, we set Zj!/ — xl®f '■— If, the identity op- 
erator. Using the defined fractional integrals (JT]) and ([2]), Riemann-Liouville fractional derivatives 
are defined as: the left Riemann-Liouville fractional derivative (RLFD) 



the right RLFD 



*ww = w ^f x j x «-*)-'*mdt = {-f x )j 1 b -"m, (4) 
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where a is the order of the derivative. The operators (H])-(EQ) are obviously linear. Below we 
present the rules of fractional integration by parts for RLFI and RLFD which are particularly 
useful for our purposes. 

Lemma 1. ([26]) Let < a < 1, p > 1, q > I, and 1/p + 1/q < 1 + a. If g G L p ([a,b]) and 
f G L a ([a,b}), then J a g(x) a I"f(x)dx = J a f(x) x I£g(x)dx. 

Lemma 2. [cf. [351] Let < a < 1 and functions f,g obey the assumptions of Lemma\l\ If 
f(a) = f(b) = (or the analogous condition for g is fulfilled), then 

b pb 

g(x) a D°f(x)dx = f(x) x D^g(x)dx. (5) 

J a 

Partial fractional integrals and derivatives are natural generalization of the corresponding oned- 
imensional fractional integrals and derivatives, being taken with respect to one or several variables. 
For (xi, . . . , x n ), (a±, . . . , a n ), where < on < 1, i — 1, . . . , n and [a\, bi] X . . . X [o n , b n ], partial 
Riemann-Liouville fractional integrals of order a k with respect to Xk are defined by 

1 f Xk 
a h Ixuf( x i>---i x n) = ^r, — r / {xk -t k ) ak ~ 1 f(x 1 ,... } x k -i,t k ,x k+1 ,...,x n )dt k , x k >a k , 
r(a fe ) J ak 

1 f bk 
x k Ib h k f( x i>---> x n) = ^7 — 7 / (** -x k ) ak ~ l f{x ll ...,x k -i,t k ,x k+ i,...,x n )dt k i x k <b k . 

- 1 - V a k) J Xk 

Partial Riemann-Liouville derivatives are defined by: 

1 d f Xk 

a k D" k f(x 1 ,...,x n ) = -— — / (Xk-tk) ak f{xi,...,x k -i,t k ,x k+ i,...,x n )dt k , 

r(l - a k ) dx k J ak 

1 d f bk 

x k D? k f(xi,...,x n ) = -— — / (tk-Xk) ak f(xi,...,x k -i,t k ,x k+ i,...,x n )dt k . 

k r(l - a k ) dx k J Xk 

3 Fractional Euler— Lagrange Equations 

Consider a physical system characterized by a set of functions 

Uj(t,xi,...,x n ), j = l,...,m, (6) 

depending on time t and the space coordinates Xi, . . -,x n . We can simplify the notation by 
interpreting ^ as a vector function u — (u%, . . . , u m ) and writing t — Xq, x — (xq,Xx, . . . ,i„), 
dx — dxodxi ■ ■ ■ dx n . Then ((5]) becomes simply u(x) and is called a vector field. Following 
[H1CEI], we introduce the operator V Q defined by V Q = ( ao D x °, ai D^, ■ ■ ■ , a „-D"^), where a = 
(ao, (Xx, ... , ct n ), < on < 1, i = 0, . . . ,n. Consider the action functional in the form 

pbo n 

J(u)= L(u,V a u)dx = £{u,V a u)dx, (7) 

J ao J O 

where R = [a\, b{\ x . . . x [a n , b n ], fl — Rx [ao, bo}. The boundary conditions are specified by asking 
that u(x) = (f(x) for x G dft, where <p : dfl — > W m is a given function. The functions L(u, V a u) 
and C(u, \7 a u) are called the Lagrangian and Lagrangian density of the field, respectively. The 
advantage of using V Q in the formulation of problem ([7]) lies in the fact that fractional derivatives 
take into account the space history and time history. Therefore, comparing with integer order 
derivatives, they are more appropriate for representation of natural phenomena. We assume that: 
(i) The set of admissible function consists of all function u : ft — > M m , such that ai D"fv,j G C(Q,, M), 
i = 0, . ..,n, j = l,...,m, and u\ dn = <p\ (ii) C G C 1 {W n x R m (" +1 );M); (iii) Xi D% gjfa^ , * = 

0, . . . , n, j = 1, . . . , m are continuous on (ao, bo) X . . . x (a„, b n ) and integrable on £1 for all admissible 
functions. Applying the principle of stationary action to we obtain the multidimensional 
fractional Euler-Lagrange equations for the field. 



Theorem 1. If function u minimizes the action functional ([7]), then u satisfies the multidimen- 
sional fractional Euler- Lagrange differential equations: 

H+t-^dl^- - 3=1 -- m - (8) 

J j_q at J-i J 

Proof. The variations for this problem are functions h : $7 — > M. m such that atD^hj £ C(f2, K), 
j = 0, . . . , n, j = 1, . . . , to, and h vanishes on the boundary: h(0) = 0, x £ dfl. The condition 
that a function u : il — >• R m be a critical point of the action »7(u) is -j-Jiu + eh)\ = for all 
variations h with h\esi — 0- Differentiating under the integral sign we obtain 
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The Fubini theorem allows us to rewrite integrals as the iterated integrals so that we can use the 
integration by parts formula (JSJ: 

J2^ D h a tyh.. )hjdx, j = l,...m 
since h\dd — 0. Therefore, 

By fundamental lemma of the calculus of variations (see Section 2.2 in [31]), u satisfies (J5J). D 

In the following we present physical examples. The main aim is to illustrate what are the 
Euler— Lagrange equations when we have a Lagrangian density depending on fractional derivatives. 

Example 1. Consider a minimizer of J(u) — ^ J J R J2i=i {{ai^xi u ) ~ f u ) dx, where f : R^R 
is a given function, in a set of functions that satisfy condition u = p on dR, where tp is a given 
function defined on the boundary dR of R = [ai, bi] x [02, 62]- By Theorem [l] a minimizer satisfies 
the following equation Yli=i Xi D" i ai D^ju = f. Observe that if on goes to f, then the operator 
ai D"*, i = 1,2, can be replaced with J-, and the operator Xi D™', i = 1,2, can be replaced with 
— jp- (see [M]). Thus, for a — s> 1 we obtain the Poisson equation that arises in the potential theory 
in electrostatic fields. 

Example 2. Consider a motion of medium whose displacement may be described by a scalar 
function u(t,x), where x — {x\,X2)- For example, this function might represent the transverse 
displacement of the membrane. Suppose that the kinetic energy T and potential energy V of 
the medium are given by: T (^) = \ J p (^) dx, V(u) = \ J k\\ Vu\\ 2 dx, where p(x) is a mass 
density and k(x) is a stiffness, both assume positive. Then, the classical action functional is 

J(u) = \ J J Q (p(§f) ~^||Vu|| 2 j dxdt. When we have the Lagrangian with the kinetic term 
depending on a fractional derivative, then the fractional action functional has the form 

J{u) = \l j (p (oA a °«) a ~ fell Vwf) dxdt. (9) 

The fractional Euler-Lagrange equation satisfied by a stationary point of © is 

p t D*° Q D*°u - V(fcVu) = 0. 



If p and k are constants, then equation tDf^Dfu — c 2 Au = where c 2 = k/p can be called 
the time-fractional wave equation. In the case when the kinetic and potential energy depend on 
fractional derivatives, the action functional for the system has the form 

J{u) = \j J [p(oD?°u) 2 - k(( ai D^u) 2 + ( a2 D%u) 2 )] dxdt. (10) 

The fractional Euler-Lagrange equation satisfied by a stationary point of (|10l) is 

PtD^ D^u - k( Xl D% ai D£u + X2 D% a2 D%u) = 0. 

If p and k are constants, then t D" ^Dfu — c 2 ( Xl D^ 1 ai D x ^u + x 2 D^ 2 a2 D x ^u) — can be called 
the space-and time-fractional wave equation. Observe that in the limit, a, ot\,(X<2. — \ 1, we obtain 
the classical wave equation 4^ — c 2 Au — with wave-speed c. 

Remark 1 . Recently, several authors have investigated the fractional diffusion/wave equation and 
its special properties. Mostly fractional wave equations are obtained by changing classical partial 
derivatives by fractional, it can be the Riemann-Liouville, Caputo or another one 6, 29, 30"! |3"TI 155] . 
Here, following |H1 HH 132] we use the different approach. We start with a variational formulation 
of a physical process in which one modifies the Lagrangian density by replacing integer order 
derivatives with fractional ones. Then the action integral in the sense of Hamilton is minimized 
and the governing equation of a physical process is obtained. This approach can be explain 
via the notion of fractional embedding, we refer the reader to [8]. The main remark is that 
equations by themselves do not have a universal significance, their form depending mostly on the 
coordinates systems being used to derive them. On the contrary the underlying first principle like 
the least-action principle carry an information which is of physical interest and not related to the 
coordinates system which is used. 

Remark 2. It should be mentioned, that since fractional operators are nonlocal, it can be extremely 
challenging to find analytical solutions to fractional problems of the calculus of variations and, in 
many cases, solutions may not exist. The question of existence of solutions to fractional variational 
problems is a complete open area of research. This needs attention. As mentioned by Young, the 
calculus of variations has born from the study of necessary optimality conditions, but any such 
theory is "naive" until the existence of minimizers is verified. The process leading to the existence 
theorems was introduced by Leonida Tonelli in 1915 by the so-called direct method. During two 
centuries, mathematicians were developing "the naive approach to the calculus of variations". 
Similar situation happens now with the fractional calculus of variations: the subject is only fifteen 
years old, and is still in the "naive period". We believe time has come to address the existence 
question, and this will be considered in a forthcoming paper. 

4 Fractional Noether-type Theorem 

In this section we prove the fractional Noether-type theorem. We start by introducing the notion 
of variational invariance. 

Definition 1. Functional © is said to be invariant under an e-parameter group of infinitesimal 

transformations u(x) = u + e£(x, u) + o(e) if 

f C{u,V a u)dx= f C{u,V a u)dx (11) 

for any fi* C il. 

The next theorem establishes a necessary condition of invariance. 
Theorem 2. If functional ([7]) is invariant in the sense of Definition^ then 



j — \ J j — I %—Q J 



Proof. Equation (fTTj) is equivalent to C(u, V Q w) = C(u + e£(x, u) + o(s), V a (u + e£(x, u) + o(e))). 
Differentiating both sides of with respect to e, then substituting e = 0, and using the definitions 
and properties of the fractional partial derivatives, we obtain equation (|12[) . □ 



The following definition is similar to 17] and useful in order to formulate the fractional Noether- 
type theorem. 



Definition 2. For < 7 < 1 and a given ordered pair of functions / and g (/,<?: R" +1 — > M) we 
introduce the following operator: V (f,g) = f • ai D1_g — g ■ XiD^. f , where xi S [aj,&i], i = 0, . . . , n. 

We now prove the fractional Noether-type theorem. 

Theorem 3. If functional |[7J) is invariant in the sense of Definition^ then 

n m / fir \ 

£i>-y§^ =° <i3) 

along any solution of ([5]), 

Proof. Using the fractional Euler-Lagrange equations ©, we have: 

J— 1 J J—l 2—0 J 



Substituting (fl4|) in the necessary condition of invariance (fT2|) we get 

o = - E E * *£■ a-^r^ + E E ^w^a 

j — \ %— J j — 1 z— J 

n / m „ „ m „„ 

= E E a D on aiDZiZj - E a! *- D 6i' « D a *u^ j 

By definition of operator I? 7 we obtain equation (fT"3")l . D 

Remark 3. In the particular case when ai goes to 1, ? = 0, . . . , n, we get from the fractional conser- 
vation law (|13[) the classical Noether conservation law without transformation of the independent 
variables (that is, conservation of current): 

along any extremal surface of 3{u) = /<-, C(u, "S7u)dx. 

Example 3. Let us consider again the Lagrangian given in Example [2] L = T — V. Then the 
Hamiltonian, H = T + V, is conserved (see [20]). However, if the kinetic and potential energy 
depend on fractional derivatives, that is the action functional has the form (fT0|) . then it can 
be showed that the Hamiltonian is not conserved (cf. [2H ES] ) - Observe that functional (fTUl) 
is invariant (in the sense of Definition [T]) under transformations u(t,x) — u(t,x) + et aa ~ 1 (x\ — 
<Xi) ai (x2 — c^)" 2 . Therefore, by Theorem[3l we obtain: 

X> Ql ( 7 r^^^o°~ 1 (^-^r i - 1 (^-a 2 r- 1 ) =o, (i6) 

i=0 \°aiJ^XiU J 

with xq = t and oq = 0, along any extremal surface of (fTU)) . Note that when a* goes to 1, i = 0, 1, 2, 
expression (fT6)) has the form (|15[) . This shows that the proved fractional Noether-type theorem 
can be applied for conservative and nonconservative physical systems. 
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